Surface acoustic wave propagation in monolayer graphene 
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Surface acoustic wave (SAW) propagation is a powerful method to investigate 2D electron systems. 
We show how SAW observables are influenced by coupling to the 2D massless Dirac electrons of 
graphene and argue that Landau oscillations can be observed as function of gate voltage for constant 
field. Contrary to other transport measurements, the zero-field SAW propagation gives the wave 
vector dependence of graphene conductivity for small wave numbers. We predict a crossover from 
Schrodinger to Dirac like behaviour as a function of gate voltage, with no attenuation in the latter 
for clean samples. 
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Sound waves in solids may be used as diagnostic tool 
to investigate the low frequency long wavelength response 
of electrons. In metals they couple via the deformation 
potential mechanism, i.e. stress dependence of kinetic 
energy. In semiconductors the piezoelectric mechanism 
which generates internal electric fields that accelerate the 
charge carriers is most effective provided the compound 
has a piezoelectrically active structure. Close to smooth 
surfaces the acoustic modes may become localized in the 
direction of the surface normal and propagate as sur- 
face acoustic waves (SAW). The most prominent ones are 
Rayleigh waves [l| where the displacement vector u(i) 
lies in the 'sagittal plane' formed by the surface normal 
vector n and the acoustic propagation vector q which is 
parallel to the surface (normal to n = z, cf. Fig. [TJ) . The 
Rayleigh SAW may be thought as consisting of longitudi- 
nal and transverse components whose amplitudes depend 
on the normal direction in a manner that the surface ful- 
fills the mechanical stress free boundary condition. The 
Rayleigh wave is therefore an elliptically polarized wave 
in the sagittal plane which decays exponentially into the 
bulk with an acoustic penetration depth that is of the 
order of the bulk wavelength. 

Because of these properties the Rayleigh SAW exhibit 
interesting effects due to the coupling to charge carriers, 
particularly in the presence of a magnetic field. Two of 
those are well known and have been explored experimen- 
tally, i) when the field is within the surface plane and 
perpendicular to the sagittal (xz) plane then due to the 
breaking of time reversal invariance the velocity of prop- 
agation of SAW with wave vectors q and -q is different. 
This non-reciprocal SAW effect has been investigated for 
metals in detail in Rcf. S, 0. ii) In semiconductor het- 
erostructures the charge carriers are confined in the lat- 
eral (xy) plane and a 2D electron gas (2DEG) forms. 
This may be nicely investigated by Rayleigh SAW which 
are localized on the surface sheet. For fields perpendicu- 
lar to the surface they probe the longitudinal wave vec- 
tor dependent conductivity of the 2DEG. Since the latter 



exhibits the typical Landau quantum oscillations under 
suitable conditions this may be directly seen in the SAW 
propagation velocity and attenuation as function of mag- 
netic field. This idea was spectacularly confirmed for the 
2 DEG in GaAs hetero junctions Extended investi- 

gations have shown that this method also allows uniquely 
to determine the wave- vector dependence of longitudinal 
conductivity [f3]. 

The 2DEG in GaAs heterostructure contains electrons 
with non-relativistic parabolic dispersion and hence a fi- 
nite band mass. Recently graphene, a single sheet hon- 
eycomb carbon layer, has been identified as a genuine 
2DEG systemQ. It contains quasi-relativistic massless 
Dirac electrons with a linear dispersion around the Dirac 
point K where the slope is given by the Fermi velocity 
vp = 10 6 m/s, obeying to 



H = V F (a x p x + (JyPy), 



(1) 



where the pseudospin variables (a) arise from the two- 
sublattice structure. The transport and thermodynamic 
properties of this system has been much discussed and 
investigated, in particular in an external field [H, 0]. 
They arc distinctly different from the non-relativistic case 
which is mostly due to two effects: i) because of the linear 
dispersion the density of states (DOS) at the Dirac point 
E = vanishes according to N{E) = \E\/D 2 (per valley 
and spin; D 2 = 2irv F /A Ci where A c is the area of the 
hexagonal cell). As a consequence the screening proper- 
ties of Dirac electrons will be different from those of the 
nonrelativistic case 1^, 11 [. ii) The Landau quantization 
in an external magnetic field leads to relativistic Landau 
levels at E a (n) = au c \/n + 1 with a = ±1, n = 0, 1, 2.. 
and the Landau scale ui c = vpy // 2e\B\. In addition there 
is a field independent level always at E* = 0. 

There has been considerable discussion whether the 
Dirac point and the linear DOS in single layer graphene 
is stable against perturbations like buckling, disorder or 
excitonic gap formation. We ignore these subtle issues 
here and refer to the literature [9f . Furthermore we note 
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FIG. 1: Schematic setup for SAW propagation on monolayer 
graphene sandwiched between piezoelectric slabs with dielec- 
tric constant e as adapted from Ref. ITM SAW transducer and 
detector are attached to the top slab with thickness qd <C 1 (q 
= SAW wave number ). The SAW displacement u(t) at fixed 
position traces an ellipse in the sagittal (xz) plane and the 
penetration depth in -z direction is ~ ^iv/q 3> d. Direction of 
field B is normal to the surface. 



that the chemical potential fi may be easily moved away 
from the Dir ac p oint by applying a gate voltage V g via 
V g ~ n ~ ^ 2 [13l| (with n the carrier density). This is a 
new aspect as compared to the degenerate nonrclativistic 
2DEG and we will employ it in our discussion of SAW 
propagation effects in graphene. Preliminary SAW ex- 
periments in zero field have been performed in Ref. [14 . 
In this setup the aim was rather to detect surface adsor- 
bates via their change of graphene electronic properties 
and ensuing SAW propagation anomalies. In the present 
work we want to look at intrinsic properties of graphene 
and investigate the field dependence, gate voltage (chem- 
ical potential) and wave vector dependence of the SAW 
propagation velocity and attenuation. A possible exper- 
imental setup which we have in mind is shown schemat- 
ically in Fig. [1] The graphene sheet is sandwiched be- 
tween two insulating piezoelectric layers with the same 



background dielectric constant e, similarly to Ref. [15. 
The SAW r.f. transducer and the detector are mounted 
on the top layer. Its thickness d is assumed to be much 
smaller than the penetration depth of the SAW along 
—z. There may also be alternatives where the coupling 
of SAW in the piezoelectric substrate to the Dirac elec- 
trons happens purely capacitive across an empty space. 

The analysis of SAW propagation in this setup has 
two distinct steps. Firstly the mechanical and electro- 
dynamic boundary value problem for the SAW has to 
be solved which gives the SAW velocity change in terms 
of the (magneto-) conductive properties of the graphene 
sheet. This problem has been analyzed and solved in 
Ref. [l2[ In a second step we calculate the magnetocon- 
ductivity of Dirac electrons as function of field and chem- 



ical potential which leads to the SAW velocity change 
and attenuation as function of field and gate voltage. In 
this part we follow the procedure described in Ref. 
Alternatively we consider the zero field case and investi- 
gate the wave vector (or frequency) dependence of SAW 
properties, which is related to the wave vector dependent 
(B = 0) polarization given in Ref. [l(J [ll|. 

If one neglects the influence of the 2D relativistic elec- 
tron gas the SAW propagation is a purely mechanical 
problem which is determined by the elastic constants. 
We don't make any specific assumptions about the bulk 
and top layer material but rather treat it as an isotropic 
piezoelectric medium characterized by longitudinal (L) 
and transverse (T) elastic constants or velocities cp or 
v = V ( c /3 / p) (P = L7T 1 ) respectively where p is the 
mass density. The stress-free boundary condition for the 
geometry in Fig. [1] is <Ji Z = (i = x,y,z). The solution 



of the boundary problem [17| leads to the velocity of the 
surface Rayleigh waves v s =£,vt- Here v s is smaller than 
the transverse bulk velocity by a factor 0.87 < £ < 0.95 
depending on the ratio cl/ct- Due to the piezoelectric 
coupling an external potential is created which couples to 
the 2DEG via the density response function. The latter 
is completely determined by conservation laws and linear 
response relations. The induced density in the 2DEG 
leads to an additional energy density which has to be 
added to the elastic energy density causing a renormal- 
ization of the elastic constants or velocity. Because of the 
dissipative component in the density response the SAW 
velocity change is accompanied by an attenuation. It 
may be seen as Landau damping due to particle-hole ex- 
citations in the graphene sheet. By the above procedure 
fl2| SAW velocity change Av s and attenuation coefficient 
K are obtained as 



Av a (q,B,fj.) iK(q,B,fi) 



<qd) 2 /2 



1 + ia xx (q,uj,B,fj,)/a ri 



(2) 

Here us = v s q is the SAW frequency and q || x the wave 
number. The penetration depths k£ , of displace- 
ment components u x , u z into the bulk along —z direction 
are of the order of the wave length A = 2n/q. The ef- 
fective coupling a(qd) of SAW to the graphene 2DEG 
is a complicated expression composed of phenomenolog- 
ical parameters which will not be given here. It de- 
pends nonmonotonically on qd and in the limit qd <C 1 
or K(l,t)Q *C 1 assumed here is simply a constant to be 
determined by experiment. 

The central quantity is the longitudinal electronic con- 
ductivity a xx (q,uj; B, fj,) where us = v s q. The normaliza- 
tion constant is a m = v s e e ff/2ir with e e // = (erj + e )/2 
for qd < 1. The typical SAW frequency of ~ 0.1 GHz 
is extremely small (~ 5 x 10~ 3 K) compared to the typ- 
ical electron hole excitation energy (given by the chem- 
ical potential \i so that one may use the static limit 
&xx{q', B, /j,) = a xx (q, us = 0; B, fx)). We will consider two 
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FIG. 2: Landau oscillations for constant field (corresponding 
to N = 1000 LL within the cutoff frequency D) as func- 
tion of the chemical potential fi tuned by the gate voltage for 
T — 0.001D, m = 0.001. Upper panel: longitudinal conduc- 
tivity. Middle/lower panel: SAW velocity change/attenuation 
for a m /a xx (fi = 0) = 10 (black dashed-dotted line), 1 (blue 
solid line) and 1/10 (red dashed line). 



different cases: i) finite field case in the presence of scat- 
terers, which allows to neglect the possible q dependence 
of the conductivity if the scattering rate exceeds VFq 
and ii) zero field ultraclean case. For the former, due 
to Landau quantization, oscillations in Av s (B,(i) will 
appear both as function of field and as function of the 
chemical potential n(V g ) or gate voltage V g . The latter 
is an attractive new feature peculiar for SAW in graphene 
but not for the non-relativistic 2DEG. One may keep the 
field constant and instead observe the Landau oscilla- 
tions as function of gate voltage. For the first case we 
need the field dependent conductivity <j xx (B,[a) which 
was derived for localized disorder (vacancies representing 
the unitary limit [3]) in the self-consistent Born approx- 
imation (SCBA) in Ref. The result is shown in the 
top panel of Fig. [2] as function of the chemical potential 
for a constant magnetic field for ui c = D/y/N + 1 with 
N = 1000. Clearly the Landau oscillations in the con- 
ductivity also influence the SAW properties. When the 
chemical potential crosses a LL, the conductivity peaks, 
and Av s develops a dip. In between LL's, the conductiv- 
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FIG. 3: Real (red dashed line) and imaginary (blue solid line) 
part of the conductivity for B = as function of chemical 
potential and wave number (inset) with ctq = e 2 v B /vF- The 
Dirac and Schrodinger like regions are separated at 2vfq = H- 



ity is suppressed, causing the enhancement of the SAW 
velocity. For large chemical potential when the graphene 
sheet becomes more metallic, both Av s and k are sup- 
pressed. 

The SAW velocity measurement in zero field presents 
a unique feature which cannot be realized by the usual 
transport measurements, which usually probe the con- 
ductivity at w > v F q, i-C. in the dynamic limit. Due 
to the finite SAW frequency the density response of the 
Dirac electrons may be probed at finite wave vector 
q = lu/v s with to <C vpq, thus probing the static prop- 
erties of the electric response. Measurement of Av s for 
various frequencies to can therefore determine the density 
response of Dirac electrons for finite wave vector to/v s . 
This has actually been proposed and carried out for the 
non-r elativistic 2DEG in Ref. 

The q (|| a;) and to dependent longitudinal conducti vity 
is calculated similarly to the polarization function (ToL 11 1 , 
and is given in the limit of to = v s q <C (p, VFq) as 



4v s e 2 

TTVF 



v s , / v F q\ , . , (v F q 



where 



/l(a .) = vi_5_ e( l 



X 2 



x arccos | — 
x 



(3) 



(4) 

e(x-i), 

(5) 



and is shown in Fig. [3l It depends only on VFq/ ' n, thus 
exhibits scaling behaviour. For small frequencies, a xx 
vanishes with to as follows from Eq. ^ via to = v s q, 
which explains the v s term on the r.h.s of Eq. The 
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FIG. 4: Schematic view of the SAW velocity change (blue 
solid line) and attenuation (red dashed line) as a function of 
the chemical potential for B — clean graphene. For small 
H, the attenuation is negligible in the Dirac regime. 



real part of the conductivity, which determines the at- 
tenuation, is suppressed by the additional v s /vp factor, 
thus the propagation velocity is more strongly influenced 
by the presence of graphene than the attenuation. 

The density response or polarizability H(q) is related 
to the longitudinal conductivity by 



%q 



(0) 



which follows from the charge continuity equation, and 
involves only the longitudinal component of the conduc- 
tivity. For q = (q x ,0), q x = q. Using the equation for 
the SAW velocity change and neglecting the attenuation 
we may easily obtain 



11(g) ~ qa n 



1 



2Av s (q) v s 



(7) 



from the measured velocity change. The density re- 
sponse for monolayer graphene has also been calculated 
in Ref. Ql], EL and SAW measurements can yield to the 
static polarizability as well. 

There are two limits to be considered depending on 
whether q is larger or smaller than twice the radius of the 
Fermi circle kp = /x/vp. For 2kp > q what we refer to as 
Schrodinger like behaviour, a xx = 8v s e 2 kp(^j — i)/irv F q, 
which holds true in a normal metal as well, and reveals 
the important kp/q dependence of the conductivity. On 
the other hand, the peculiarity of graphene is the tun- 
ability of its Fermi energy with a gate voltage, which 
provides access to the 2kp < q region, termed Dirac 
like regime, in which case <tq = —iev s /vp exactly at 
the neutrality point (fi = 0). Thus, there is practically 
no attenuation in this regime, only velocity renormaliza- 
tion. Typical SAW wavelengths range from 2-150 ^tmfioj]. 
which are translated to vpq ~ 2 — 200 K using Vp for 
graphene. Graphene samples being ballistic over the 
micron-submicron scale imply a scattering rate around 
50 K, thus, the above Dirac to Schrodinger like crossover 
can in principle be observed in clean samples with a tiny 
residual scattering rate at the neutrality point. The re- 
sulting behaviour of SAW observables are shown in Fig. 



IH which resembles closely to that of the bare conduc- 
tivity. We also speculate that SAW at // = can iden- 
tify the role inhomogeneities extending over several lat- 
tice constant (e.g. ripples, puddles) in graphene. When 
the SAW wavelength becomes comparable to their spatial 
extension, all waves passing through are scattered, and 
the sample is not transparent any more but acoustically 
opalescent. 

In summary we have studied SAW propagation in 
monolayer graphene. By exploiting the chemical poten- 
tial dependence of its conductivity, Landau oscillations 
are predicted in both the SAW velocity change and atten- 
uation by varying the gate voltage, a feature which is only 
accessible due to the quasi-relativistic excitation spec- 
trum in graphene, clearly absent in a normal 2DEG. SAW 
measurements in graphene can provide us with a deeper 
insight into the unconventional nature of the quantum 
Hall effect. Without magnetic field, the SAW frequency 
changes measure directly the wavevector dependence of 
the longitudinal conductivity, which can reveal a Dirac 
to Schrodinger like crossover, with no attenuation in the 
Dirac regime. In addition, SAW propagation in pris- 
tine graphene can yield significant information about the 
dominant scattering mechanisms close to the neutrality 
point through cr xx (q x ). 
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